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Abstract 

We prove fluctuation bounds for the particle current in totally asym- 
metric zero range processes in one dimension with nondecreasing, concave 
jump rates whose slope decays exponentially. Fluctuations in the charac- 
teristic directions have order of mag nitudet^/^. This is in agreement with 
the expectation that these systems lie in the same KPZ universality class 
as the asymmetric simple exclusion process. The result is via a robust 
argument formulated for a broad class of deposition-type processes. Be- 
sides this class of zero range processes, hypotheses of this argument have 
also been verified in the authors' earlier papers for the asymmetric simple 
exclusion and the constant rate zero range processes, and are currently 
under development for a bricklayers process with exponentially increasing 
jump rates. 
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1 Introduction 

This paper studies anomalous current fluctuations of attractive interacting sys- 
tems in one dimension with one conserved quantity. The family of models 
considered includes the asymmetric exclusion, the zero range, misanthrope-type 
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and many other processes. In the asymmetric case (to be specified later) the 
Eulerian scaHng of such a system leads to a (deterministic) hyperbolic conser- 
vation law with a hydrodynamic flux function Ti.{g). The characteristics of the 
conservation law is of particular importance both for the PDE itself and for the 
underlying stochastic system. Recently, the current fluctuations through the 
characteristic lines drew much attention. The behavior of these fluctuations is 
fundamentally determined by the form of Ti.. Rigorous results exist for examples 
that fall in two categories. 

Order t^/* fluctuations. When Ti. is linear the fluctuations are of order t^^^ 
and converge to Gaussian processes related to fractional Brownian motion. This 
has been proved for independent particles [T71 [251 132] and the random average 
process [8j[19]. 

Order t^/'^ fluctuations. When H"{q) ^ the fluctuations are of order 
t^/'^ and converge to distributions and processes related to the Tracy- Widom 
distributions from random matrix theory. The most-studied examples are the 
totally asymmetric simple exclusion process (TASEP), the polynuclear growth 
model (PNG) and the Hammersley process. Two types of mathematical work 
should be distinguished. 

(a) Exact limit distributions have been derived with techniques of asymptotic 
analysis applied to determinantal representations of the probabilities of interest. 
Most of this work has dealt with particular deterministic initial conditions, and 
the stationary situation has been less studied. The seminal results appeared in 
[3] for the last-passage version of the Hammersley process and in [22] for the 
last-passage model associated with TASEP. Current fluctuations for stationary 
TASEP were analyzed in [20]. Here is a selection of further results in this 
direction: [H [TH HT] US] [55] . Recently, the asymmetric simple exclusion (ASEP) 
also got within reach of these techniques [34] . 

(b) Probabilistic approaches exist to prove fluctuation bounds of the correct 
order. The seminal work [TS] was on the last-passage version of the Hammersley 
process, and then the approach was adapted to the last-passage model associated 
with TASEP ^ . The next step was the development of a proof that works for 
particle systems: the ASEP was treated in [l3] and the totally asymmetric zero 
range process (TAZRP) with constant jump rate in [7]. The ASEP work [l^ was 
the first to prove t^^^ order of fluctuations for a process where particle motion 
is not restricted to totally asymmetric. Resolvent methods were also applied in 
[30] [29] to extend the results from nearest neighbor ASEP to exclusion processes 
with non nearest neighbor jumps. 

The present paper takes a further step toward universality of the t^^^ order 
of fluctuations in the case 'H"{g) ^ 0. We rewrite our earlier proof for ASEP 
and constant rate TAZRP in a fairly general way, extract and formulate in 
general terms a particular feature of these two models that made our proof work. 
For reasons to be explained later we call this feature microscopic concavity. 
With this notion in hand we extend the t^/'^ scaling result for a class of totally 
asymmetric zero range processes (with non constant rates). We remark at this 
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point that jump rates of this example have a much richer behavior than the 
constant rates of those featured in anomalous scaling proofs so far. Further 
generalizations now only require the verification of microscopic concavity. 

The hypothesis of microscopic concavity consists of control of second class 
particles that is a microscopic counterpart of the macroscopic effect that con- 
cavity of Ti has on characteristics. We make this technically precise in Section 
12.61 Once the microscopic concavity assumption is made the proof works for the 
entire class of processes. This then is the sense in which we take a step toward 
universality. As a bi-product, we also obtain superdiffusivity of the second class 
particle in the stationary process. 

Earlier proofs of i^^'^ fluctuations have been quite rigid in the sense that 
they work only for particular cases of the models where special combinatorial 
properties emerge as if through some fortuitous coincidences. There is basically 
no room for perturbing the rules of the process. By contrast, the proof given 
in the present paper works for the whole class of processes. The hypothesis of 
microscopic concavity that is required is certainly nontrivial. But it does not 
seem to rigidly exclude all but a handful of the processes in the broad class. The 
estimates that it requires can probably be proved in different ways for different 
subclasses of the processes. And the proof itself may evolve further and weaken 
the hypothesis required. 

To summarize, we are currently able to verify the required hypothesis of 
microscopic concavity for the following three subclasses of processes. 

(i) The asymmetric simple exclusion process (ASEP). Full details of this case 
are reported elsewhere [12] and we give a brief informal description in Section 
12.8.11 This proof is somewhat simpler than the earlier one given in [13] . 

(ii) Totally asymmetric zero range processes (TAZRP) with a concave jump 
rate function whose slope decreases geometrically, and may be eventually con- 
stant. This example has been out of reach for existing methods, so it is com- 
pletely new in this context. It is developed fully in the present paper. As a 
special case, the result of [7] for the constant rate TAZRP is also recovered. 

(iii) The totally asymmetric bricklayers process with convex, exponential 
jump rate. This system satisfies the analogous microscopic convexity. Due to 
the fast growth of the jump rate function this example needs more preliminary 
work than was sensible to include in the present paper, and so the result will 
be published separately in the future. 

We expect that a broader class of totally asymmetric concave zero range 
processes should be amenable to further progress because a key part of the 
hypothesis can be verified, and only a certain tail estimate is missing. We 
explain this in Section r2. 8. 21 

Interacting particle systems can naturally be given a surface growth repre- 
sentation where integrated particle current becomes the height of a surface and 
particle occupations become (negative) discrete gradients of this surface. We 
found this picture extremely helpful in visualizing currents and couplings, hence 
this is the way we introduce and handle the processes. 

This paper has two parts. In the main part we prove the general fluctuation 
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bound under the assumptions needed for membership in the class of processes 
and the assumption of microscopic concavity. The remainder of the paper shows 
that the assumptions required by the general result are satisfied by a class of 
zero range processes. Here is a section by section outline. 

In Section [2] we define the general family of processes under consideration, 
describe the microscopic concavity property and other assumptions used, and 
state the general results. Partly as corollaries to the fluctuation bound along 
the characteristic we obtain a law of large numbers for a second class parti- 
cle and limits that show how fluctuations in non-characteristic directions on 
the diffusive scale come directly from fluctuations of the initial state. Section 
12.81 describes two examples. Section 12.8.11 gives a brief description of how the 
asymmetric simple exclusion process (ASEP) satisfies the assumptions of our 
general theorem. (Full details for this example are reported in [12].) Section 
12.8.21 describes a class of totally asymmetric zero range processes with concave 
jump rates that increase with exponentially decaying slope. 

The general theorem is proved in two parts: the upper bound in Section [3] 
and the lower bound in Section 01 Section [5] proves a strong law for the second 
class particle, partly as a corollary of the main fluctuation bounds. We then 
return to the zero range example and give a complete proof for this class of 
processes in Section [51 

A three-part Appendix contains auxiliary computations for the stationary 
distribution and hydrodynamic flux function. In particular, if the jump rate 
function of a zero range process is concave and not linear then the hydrodynamic 
flux Ti. satisfies ^{"{q) < for all densities < g < oo. 

Notation 

We summarize here some notation for easy reference. Z+ = {0, 1, 2, . . . }, K+ = 
[0, oo). Centering a random variable is denoted hy X = X — EX. Constants 
C., a, do not depend on time, but may depend on the density parameter g and 
their values can change from line to line. The numbering of these constants is of 
no particular significance and is meant only to facilitate following the arguments. 

2 Definitions and results 

We define the class of processes studied in this paper, give a list of examples, and 
discuss some of basic properties. Then come the hypotheses and main results 
of this paper, followed by two examples of subclasses of processes for which the 
hypotheses can be verified. 

2.1 A family of deposition processes 

The family of processes we consider is the one described in [TT] , and we repeat 
the definition here. We start with the interface growth picture, but we end up 
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using the height and partiele languages interchangeably. For extended-integer- 
valued boundaries — oo < a;™™ < and 1 < uj'^^^ < oo define the single-site 
state space 

I:={zeZ : -l<z< w'"'*^ + 1} 
and the increment configuration space 

fl:={uj_ = {uJi)ie-L ■■ oJiGl} = I^. 

At times it will be convenient to have notation for the increment configuration 
S^gCI with exactly one nonzero entry equal to 1: 



(2.1) {S,h 



1, fori=j, 
0, for i j. 



For each pair of neighboring sites i and i + 1 of Z imagine a column of bricks 
over the interval (i, i + 1). The height hi of this column is integer-valued. The 
components of a configuration a; e O are the negative discrete gradients of the 
heights: tOi = hi-\ — hi G I. 

The evolution is described by jump processes whose rates p and q are non- 
negative functions on / x /. Two types of moves are possible. A brick can be 
deposited: 

/oo^ (wi>Wi-n) — > {uji ~ 1, uji+i + 1)\ , . 

(2.2) > with rate ^(Wj, Wj+i), 
hi — > hi + 1 J 

or removed: 

{(Ji,U}^+l) > (UJ, + 1, U)i+1 - 1)\ ... , , , 

(2.3) , , ^ > with rate ^'(wj, Wi+i). 
hi — > hi — 1 I 

Conditionally on the present state, these moves happen independently at all 
sites i. We can summarize this information in the formal infinitesimal generator 
L of the process w(.): 



(i<^)(w) = ^p(wi, oji+i) ■ [(p{. . . , - 1, oji+i + 1, ...)- (piw)] 
+ ^ q{uJi, uii+i) ■ [(p{. . . , + 1, uii+i - 1, . . . ) - <^(w)] . 



(2.4) 



L acts on bounded cylinder functions (p : Q M. (this means that tp depends 
only on finitely many w^-values). 

Thus we have a Markov process {oj{t) : t G R+j of an evolving increment 
configuration and a Markov process {h{t) : t E K+} of an evolving height 
configuration. The initial increments w(0) specify the initial height h{0) up to 
a vertical translation. We shall always normalize the height process so that 
/io(0) = 0. 
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In the particle picture the variable uii {t) represents the number of particles 
at site i at time t. Step \2.2\ represents a rightward jump of a particle over 
the edge (i, i + 1), while step (|2.3p represents a leftward jump. (If negative 
w-values are permitted, one needs to consider particles and antiparticles, with 
antiparticles jumping in the opposite direction.) Figure [T] shows a configuration 
and a possible step with both walls and particles. It is in the particle guise that 
many of these processes appear in the literature: simple exclusion processes, 
zero range processes and misanthrope processes are examples included in the 
class studied in this paper. 



3=2 



Figure 1: The wall and the particles with a possible step 



It will be useful to see that 

hi{t) ~ hi{t)~ho{0) — the net number of particles that have passed, 
(2.5) from left to right, the straight-line space-time path 

that connects (1/2, 0) to {i + 1/2, t). 

In particular, height increment hi{t) — hi{0) is the cumulative net particle current 
across the edge {i, i + 1) during time (0, t]. 

We impose the following four assumptions ()2.6p - (|2.9|) on the rates. 

• The rates p, q : I x I ^ M+ must satisfy 

(2.6) p(c.™, • ) ^ P( • , ^ • ) ^ 9( • , ^'"'") = 

whenever either w""" or uj"-^'^^ is finite. Either both p and q are strictly 
positive in all other cases, or one of them is identically zero. The process 
is called totally asymmetric if either g = or p = 0. 

• The dynamics has a smoothing effect when we assume the following mono- 
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tonicity: 



(2.7) 



p{z + 1, y) >p{z, y), 
q{z + l,y)< q{z, y), 



p{y, z + 1) <p{y, z) 
q{y, z + l)> q{y, z) 



for y^ z, z + 1 G /. Under this property the higher the neighbors of a 
column, the faster it grows and the longer it waits for a brick removal, on 
average. This is the notion of attractivity. 

• The next two assmrrptions guarantee the existence of translation-invariant 
product-form stationary measures. (Similar assumptions were employed 
by Cocozza-Thivent [T6'.) 

— For any x, y, z d I 

p{x, y)+p(y, z) +p{z, x) 
(2.8) +q{x, y) + q{y, z) + q{z, x)^p{x, z)+p{z, y)+p{y, x) 

+ q{x, z) + q{z, y) + q{y, x). 

— There are symmetric functions Sp and Sq on I x /, and a function 
f on I such that /(w™'") = whenever w™" is finite, f{z) > for 
z > w™'", and for any y, z £ I, 



(Interpret Sp{y, z) = Sq{y, z) — if y or z > ui"'"^^.) Condition (|2.7p 
implies that / is nondecreasing on /. 

An attempt at covering this broad class of processes raises the uncomfortable 
point that there is no unified existence proof for this entire class. Different 
constructions in the literature place various boundedness or growth conditions 
on p and q and the space J, and result in various degrees of regularity for the 
semigroup. (Among key references are Liggett 's monograph [2^, and articles 
[T], [S] and [26j.) These existence matters are beyond the scope of this paper. 
Yet we wish to give a general proof for fluctuations that in principal works for all 
processes in the family, subject to the more serious assumptions we explain in 
Section [^31 To avoid extraneous technical issues we make the following blanket 
assumptions on the rates p and q to be considered. 

• We assume that the increment process Lu{t) , and the corresponding height 
process h{t) with normalization /io(0) = 0, that obey Poisson rates p and 
q as described by (|2.2p and (|2.3p . can be constructed with cadlag paths 
in a subspace fl of tempered increment configurations (i.e. configurations 
that obey some restrictive growth conditions). 

• The subspace ft is of full measure under the invariant distributions /i^ 
defined in Section [2^ 



(2.9) 



p{y, z) = Sp(y, z + l)f[y) 
and q{y, z) = Sq{y + 1, z)f{z). 
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• It is also possible to construct jointly ^everal versions of the process with 
initial configurations from the space O and with joint evolution obeying 
basic coupling (described in Section . 

• Rates p and q have all moments under the invariant distributions /i^. In 
fact arguments like Lemma lC.2l of the Appendix provide this when / does 
not grow faster than exponential on Z+ and does not decrease faster to 
zero than exponential on . 

The reader will see that our proofs in Sections [31 [H [5] and [H do not make 
any analytic demands on the semigroup and its relation to the generator. We 
only use couplings, counting of particle currents and simple Poisson bounds. 

Two identities from article PJj play a key role in this paper, given as (|2.19p 
and (|2.20p in Section [^75l These identities hold for all processes in the family 
under study. The proofs given in [11] use generator calculations which may not 
be justified for all these processes. However, these identities can also be proved 
by counting particles and taking limits of finite- volume processes (pj] contains 
an example) . Such a proof should be available with any reasonable construction 
of a process. Hence we shall not hesitate to use the results of ^11] . 

2.2 Examples 

To give concrete meaning to the general formulation of the previous section we 
describe some basic examples. The type of state space / distinguishes three cases 
that we call generalized exclusion, misanthrope and bricklayers processes. In all 
cases there are two parameters < p, q < I such that p + q = 1. Asymmetric 
processes have p ^ q. These are the processes for which our results are relevant. 

1. Generalized exclusion processes. These are the cases where both w™™ 
and Lj™^'' are finite. 

• The asymmetric simple exclusion process (ASEP) introduced 
by F. Spitzer [33] is defined by w""" = 0, uj'^'^'' = 1, f{z) = l{z = 1}, 
Sp{y, z) = p ■ l{y = z = 1} and s,(y, z) = q ■ l{y = z = 1}. This 
produces the familiar rates 

P{y, z) = p ■ l{y = I, z = 0} and q{y,z) = q-l{y = 0,z = l}. 

Here LUi G {0, 1} is the occupation number for site i, p{uji, w^+i) is 
the rate for a particle to jump from site i to i + I, and (/(wi, uji+i) is 
the rate for a particle to jump from site i + 1 to i. These rates have 
values p and q, respectively, whenever there is a particle to perform 
the above jumps, and there is no particle on the terminal site of the 
jumps. Conditions (j2.7p and (|2.8p are also satisfied by these rates. 

• Particle-antiparticle exclusion process. Let uj™™ =1, uj™'^'^ = 
1. Take /(—I) = 0, /(O) = c (creation), /(I) = a (annihilation) 
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where c and a are positive rates with c < a/2, 

sp(0, 1) = 0) = p, sp(0, 0) = 1) = |, 

5,(0, 1)= 5,(1,0) 5,(0, 0) = g, 5,(1, 1) = I 

and 5p, 5, zero in all other cases. These result in rates 

p(0, 0) = PC p(0, -1) = p{l, 0) = ^, p(l, -1) = pa, 

g(0, 0)=qc, 9(-l, 0) =<z(0, 1) - ^, g(-l, l) = qa 

and zero in all other cases. If uji is the number of particles at site 
i, with = — 1 meaning the presence of an antiparticle, then this 
model describes an asymmetric exclusion process of particles and 
antiparticles with annihilation and particle-antiparticle pair creation. 
These rates also satisfy our conditions. 

One can imagine other generalizations with bounded numbers of particles 
and/or antiparticles per site. 

2. Generalized misanthrope processes have lo™™ > —oo, uj™'^^ = oo. 

• Zero range process. Take u;™™ = 0, u;"^^^ = cx), an arbitrary 
nondecreasing function / : Z+ M+ such that /(O) — 0, 

Sp{y, z)=p and s,(y, z) = q, 

P{y, z) = pf{y) and q{y, z) = g/(z). 

Again, loi represents the number of particles at site i. Depending on 
this number, a particle jumps from i to the right with ratep/(wi), and 
to the left with rate qf{uji). These rates trivially satisfy conditions 

(EH) and am . 

3. General deposition processes have a;™'" = — oo and cj™^'' — oo. The 
height differences between adjacent columns can be arbitrary integers. 
Antiparticles are needed for a particle representation of the process. 

• Bricklayers process. Let / : Z ^ M+ be non-decreasing and 
satisfy 

/(z) • /(I - z) = 1 for all z e Z. 

The values of / for positive z's thus determine the values for non- 
positive z's. Let 

which results in 

P{y, z) ^pf{y)+pf{-z) and q{y,z)^qf{-y)+qf{z). 
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The following picture motivates the name bricklayers process. At 
each site i stands a bricklayer who lays a brick on the column to his 
left at rate pf{—uji) and on the column to his right at rate pf(uji). 
Each bricklayer also removes a brick from his left at rate qf{uji) and 
from his right at rate qf{—uJi). Conditions (|2.7p and (|2.8p hold for 
the rates. 

These were examples for which our theorem holds, provided the hypotheses 
on microscopic concavity to be described below can be verified. 

2.3 Basic coupling 

In basic coupling the joint evolution of n processes m = l,...,n, is 

defined in such a manner that the processes "jump together as much as possi- 
ble." The joint rates are determined as follows, given the current configurations 
iv}, Lv^, . . . , € 17. Consider a step of type (|2.2p over the edge (i, i + 1). Let 
m i{m) be a permutation that orders the rates of the individual processes 
for this move: 

r(m)=p(w/ ' , LU^^-^') < p[uj^' I ) = r(m + l), l<m<n. 

Set also the dummy value r(0) = 0. Now the rule is that independently for each 
m — 1, . . . , 71, at rate r(m) — r(m — 1), precisely processes ^ . . . ^ 

^^(n) execute move (|2.2p . and processes u/^^\ ^^^"^K ■ ■ ■ , w^*^™"^^ do not. The 
combined effect of these joint rates creates the correct marginal rates, that is, 
process tj^^™) executes this move with rate r(m). 

Notice also that, due to (|2.7p . a jump of u/" without tJ' can only occur if 
p{LJi, ujt+i) < piuJi, i^f+i) which imphes Luf > uj\ or ljj°;^t^ < uji_^^i. The result 
of this step (12. 2p then cannot increase the number of discrepancies between the 
two processes, hence the name attractivity for (12. 7p . In particular, a sitewise 
ordering ojf < lo\ Vi G Z is preserved by the basic coupling. 

One can check that moves of type (|2.3p with rates q obey the same attrac- 
tivity property. 

The differences between two processes are called second class particles. Their 
number is nonincreasing. In particular, if cj" > for each i E then the 
second class particles are conserved. In view of (|2.5I) . in this case the net 
number of second class particles that pass from left to right across the straight- 
line space-time path from (1/2, 0) to (i + 1/2, t) equals the growth difference 

(2.10) (hnt) - /iS(O)) - (h'iit) - h^O)) = h'^it) - hUt) 

between the two processes w''(-) and 'rii'i )- 

A special case that is of key importance to us is the situation where only 
one second class particle is present between two processes. 
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2.4 Translation invariant stationary product distributions 

The results of this paper concern stationary processes with particular product- 
form marginal distributions that we define in this section. For many cases it 
has been proved that these measures are the only extremal translation-invariant 
stationary distributions. Following some ideas in Cocozza-Thivent [16] . we first 
consider the nondecreasing function / whose existence was assumed in (j2.9[) . 
For I 3 z > define 

z 

/(z)!:=n/(2/), 



while for / 9 z < let 



n fiv) 

y=z+l 



and then /(O)! : = 1. This definition satisfies /(z)! • f{z + 1) = f{z + 1)! for aU 
z G /. Let 



log ( liminf (/(z)!)"'^^ ) = lim log(/(z)), if a; 

V z — >oo / z — *oo 

oo, else 



max 



and 



log ( limsup(/(-z)!)"^/^ ) = lim log(/(-z)), if tj""" = -oo 

V z->oo J 

— oo, else. 

By monotonicity of /, we have > 9. The case 9 ~ 9_ would imply that 
^min ^ ^max ^ ^^^^ j -g coustaut. Noticc that ((27)) and ((2^ imply 

that Sp is non-increasing in its variables, but p is non-decreasing in its first 
variable. Hence a constant / results in an Sp that does not depend on its first 
variable. But then by its symmetric property it does not depend on its second 
variable either, and we conclude that a constant / implies constant rates p (and, 
similarly, q). We exclude this uninteresting case by postulating 

(2.11) Assume / to be such that 9< 9. 
For 9 £ [9_, 9) define the state sum 

(2.12) zw..i:^<». 

Let the product-distribution on Q — have marginals 

(2.13) /(z) = / : c.. = z} : = ^ . (z G /). 
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Assumptions dMl), (EZl), (EH), imply that for Be (9,0) the product 

distribution fi^ is stationary for the process generated by (|2.4[) (see [11]). For 
some calculations in the Appendix it will be convenient to note that the family 
{ijP} can be obtained by exponentially weighting a probability measure /i^" for 
a fixed value 60 G (0, 9). 

P^, E^, Var^, Cov^ will refer to laws of a process evolving in this stationary 
distribution. In the appendix we show that the density 

g{9) ■. = E'{uj) 

is a strictly increasing, infinitely differentiable function of the parameter 9 that 
maps the interval {9_, 9) onto the interval (w""", w™'*'^). (The following point 
should cause no confusion: the single-site state space / consists of the integers 
between w""" and o;'"^'^, including endpoints if finite, but for density values the 
interval (w™'", uj^^^) is an interval of real numbers.) For most cases we shall 
use the density g, rather than 0, for parameterizing the stationary distributions. 
Accordingly, /i^, P^, E^, Var^, Cov^ will refer to laws of a density stationary 
process. 

2.5 Hydrodynamics and some exact identities 

The hydrodynamic flux is defined as 

(2.14) Hig) : = Ee{p{u;o, c^i) - q{uJo, ^1)). 

'H{g) is the expected net rate at which a given column grows, or at which 
particles pass any fixed lattice edge from left to right in a stationary density-p 
process. We show smoothness of H in Section[C]of the Appendix. It is expected, 
and in many instances proved, that asymmetric members of our class satisfy the 
conservation law 

drgiT, X) + dxn{g{T, X)) = 

in the Eulerian-scaled time and space variables T and X, see e.g. Rezakhan- 
lou [31j or Bahadoran, Guiol, Ravishankar and Saada [2]. The characteristic 
speed is the velocity with which small perturbations of the solution of this PDE 
propagate, and is given by 

(2.15) V^-.^n'ig). 

A particular expectation we shall need several times is 

(2.16) Wihiit)) ^n{g)t- gi, t>0,ieZ. 

For i = this follows from (|2.5p . and in general from the i = case together 
with ujj{t) = hj-i{t) — hj(t). 

When a stationary process is perturbed by adding a second class particle at 
the origin at time zero, we obtain two processes, lu^{-) and uj_{-). It is not a 
priori clear what the initial joint distribution of the occupation variables uj^{0), 
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a;o(0) should be. For ASEP there is no ambiguity due to the simplicity of the 
single-site state space: the only way to have a discrepancy is to set W(7(0) = 0, 
a;o(0) = 1. A useful generalization of this distribution to the broader class of 
processes involves the following family of probability measures on / introduced 
in [nj: 

^ °' z=y+l 

An empty sum is zero by convention and so if lo™^^ < oo, = 0. 

Consequently there is room for an additional particle under the p,^ distribution, 
in the sense that if w ~ then also u + 1 e I. 

To our knowledge these distributions Jl^ do not possess any invariance prop- 
erties. Their virtue is that they make identities (|2.19p and (|2.20p below true. We 
show in Section |B] of the Appendix that both /i^ and /I^ are stochastically mono- 
tone in the density g. (There is, however, no stochastic domination between fi^ 
and /I^ in general.) 

Denote by E the expectation w.r.t. the evolution of a pair (a;^(-), w(-)) 
started with initial data (recall (|2.ip ) 

(2.18) w-(0)=w(0)-^o~ ((g)Ai^) (^/i^ 

and evolving under the basic coupling. This pair will always have a single 
second class particle whose position is denoted by Q{t). In other words, iil~{t) = 
^{t) ~ ^Q(t)- Corollaries 2.4 and 2.5 of ll] state that 

(2.19) Yar^ihiit)) = Var^(w) • E|Q(t) - i\ 
and 

(2.20) F,{Q{t)) = ■ t 

for any z e Z and t > 0. Note in particular that in (|2.19p the variances are taken 
in a stationary process, while the expectation of Q{t) is taken in the coupling 
with initial distribution (|2.18p . 

2.6 Microscopic concavity 

From now on fix the jump rates p, q : I x I that define the process in 

question, assumed to satisfy all the assumptions discussed thus far. The t^/"^ 
current or height fluctuations are expected when the hydrodynamic flux Ti.{g) 
is strictly concave or convex. In this paper we discuss only the concave case. 
Concavity implies that the characteristic speed ='H'{q) is a nonincreasing 
function of density g: 

(2.21) \<g =^V^>V''. 
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The microscopic counterpart of a characteristic is the motion of a second 
class particle. Our key assumption that we term microscopic concavity is that 
the ordering (|2.2ip can also be realized at the particle level as an ordering be- 
tween two second class particles introduced into two processes at densities A and 
Q. Since this is now a probabilistic notion, there are several possible formula- 
tions, ranging from almost sure {Q^{t) > Q^{t) in a coupling) to distributional 
formulations. Assumption 12.11 below gives the precise technical form in which 
this paper utilizes this notion of microscopic concavity. It stipulates that the 
ordering of second class particles is achieved by processes that evolve on the 
labels of auxiliary second class particles, and also requires some control of the 
tails of these random labels. 

We do not imagine that this precise formulation will be the right one for 
all processes. We take it as a starting point and future work may lead to 
alternative formulations. Assumption 12. ll has the virtue that its requirements 
can be verified for some interesting processes. 

Let \ < g he two densities. Proposition IB. 41 in the Appendix gives the 
stochastic domination < . Define /I^ -I- 1 as the measure that gives weight 
— 1) to an integer z such that w™'" < z < uj^'^^ -f 1. Let fi^'^ be a coupling 
measure with marginals fl^ and /i^ -|- 1 and with the property 

(2.22) /^^'''{(y, z) : w™" -l<y<z< uj"'^'' + 1} = 1. 

Let also /i^'^ be a coupling measure of site-marginals /i^ and of the invariant 
distributions, with 

(2.23) M^'^'Ky, z) : w'"'" -l<y<z< uj"'^'' + 1} = 1, 

this is possible by Corollary IB. 31 of the Appendix. Note the distinction that 
under 'fl^-^ the second coordinate is strictly above the first. 

To have notation for inhomogeneous product measures on let A = (Ai)igz 
and g = {giji^z denote sequences of density values, with A^ and gi assigned to 
site i. The product distribution with marginals fl^^'Bo at the origin and 
at other sites is denoted by 

(2.24) ^^'^ := ((g) (g)/!^"'^". 

Measure /I-'- gives probability one to the event 

{(!7(0), ct;(0)) : 7/o(0) < wo(0), and 77,(0) < c^,(0) for ^ i G Z}. 

The initial configuration (?7(0), tii(O)) will always be assumed a member of this 

set, and the pair process {"^{1), ui{t)) evolves in basic coupling. In general 'jl-- 
is not stationary for this joint evolution. 

The discrepancies between these two processes are called the uj — rj (second 
class) particles. The number of such particles at site i at time t is uji(t) — rii(t). 
In the basic coupling the uj — rj particles are conserved, in the sense that none 
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are created or annihilated. We label the u — rj particles with integers, and let 
Xm{t) denote the position of particle m at time t. The initial labeling is chosen 
to satisfy 

■ • • < ^-i(0) < XoiO) = < Xi(0) < • • • . 

We can specify that Xq (0) = because under /t- - there is an w — r; particle at 
site with probability 1. During the evolution we keep the positions Xi{t) of 
the Lo — rj particles ordered. To achieve this we stipulate that 

whenever an oj — r/ particle jumps from a site, 
(2.25) if the jump is to the right the highest label moves, 

and if the jump is to the left the lowest label moves. 

Here is the precise form of microscopic concavity for this paper. The assump- 
tion states that a certain joint construction of processes (that is, a coupling) 
can be performed for a range of densities in a neighborhood of a fixed density 
Q. Recall p.ip for the definition of the configuration 5_. 

Assumption 2.1. Given a density g £ {uj™^"', a;™"^), there exists 70 > such 
that the following holds. For any A and g such that g — Jo < Xi < Qi < g + jo 
for all i £ Tj, a joint process {r]{t), uAt), y{t), z{t))t>o can be constructed with 
the following properties. 

• Initially (?/(0), lj{0)) is ^.-^--distributed and the joint process {ri{-), lo{-)) 
evolves in basic coupling. 

• Processes y{-) and z{-) are integer-valued. Initially y(0) = z(0) = 0. With 
probability one 

(2.26) y{t) < z{t) for all t > 0. 

• Define the processes 

(2.27) u;-it):^ui{t)~Sx^^^^^,^ and (t) : = vit) + ^x^,,,it)- 
Then both pairs (rj, 77+) and (cj~, u) evolve marginally in basic coupling. 

• For each 7 G (0, 70) and large enough t > there exists a probability 
distribution i/^''^{t) on Z+ satisfying the tail bound 

(2.28) iy'-''it){y ■■y>yo}< Ct^-'i^^-^y^^ 

for some fixed constants 3/2 < < 3 and C < 00, and such that if 
g — "f<Xi<gi<g + "f for all « G Z, then we have the stochastic bounds 

(2.29) y{t) < i/'?'''(t) and z{t) > ~v<'-'^{t). 
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Let us clarify some of the details in this assumption. 

Equation l|2.27|) says that Q^{t) :— Xz{t){t) is the single second class par- 
ticle between i] and rj^, while Q{t) : — Xy(t){t) is the one between and u. 
The first three bullets say that it is possible to construct jointly four processes 
(77, 77"*", w", Lu) with the specified initial conditions and so that each pair (77, uf), 
(ry, 77+) and w) has the desired marginal distribution, and most importantly 
so that 

(2.30) Q'^it) = > Xy^,)it) = Q{t). 

This is a consequence of (|2.26p because the u — ij particles Xi{t) stay ordered. 

The tail bound (|2.28p is formulated in this somewhat complicated fashion 
because this appears to be the weakest form our present proof allows. In our 
currently available examples i^^'^(t) is actually a fixed geometric distribution. 
However, we expect that other examples will require more complicated bounds 
and so including this generality is sensible. 

The assumptions made imply 77(i) < io{t) a.s., and by (I2.27|) 

Vit) < V^it) < ^(i) and ^(t) < uTit) < uj{t) a.s. 

In our actual constructions of the processes rj, rj^ , cj^, for ASEP (Section 
12.8.11 and [12]), for a class of totally asymmetric zero range processes (Section 
[S]) and for the totally asymmetric bricklayers process with exponential rates 
(future work) it turns out that the triples (77, 77+, lo) and (77, cj", w) evolve also 
in basic coupling, but the full joint evolution (77, 77+, lu^, lo) does not. 

As already explained, the microscopic concavity idea is contained in in- 
equality (|2.26p . There is also a sense in which the tail bounds (|2.29|) relate to 
concavity of the fiux. Consider the situation Xi = X < g = gi. We would expect 
the Lu ^ 1] particle Xq{-) to have average and long-term velocity 

n{g)-H{X) 

R[X, g) = , 

g-X 

the Rankine-Hugoniot or shock speed. By concavity Ti-'ig) — < R{X, g) < 
— H'{X). A strict microscopic counterpart would he y{t) < < z{t). But 

this condition is overly restrictive. The only cases we know to satisfy it are the 

totally asymmetric simple exclusion process and the totally asymmetric zero 

range process with constant rate. The distributional bounds (|2.29p are natural 

relaxations oi y{t) < < z{t). 

By the same token, perhaps the way to covering more examples with our 

approach involves a similar distributional weakening of (|2.26p . but this seems 

less straightforward. 

2.7 Results 

We need a few more assumptions and then we can state the main result. Con- 
stants C., a, will not depend on time, but might depend on the density pa- 
rameter g, and their values can change from line to line. We are now working 
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with a fixed member of the class of processes described in Section 12.11 with rate 
functions p, g : J x / — K"*". Recall that TL is the hydrodynamic flux defined 
in (|2.14p . In the Appendix we show Ti. is infinitely differentiable under the 
restrictions on the rates placed in Section [TT] 

Assumption 2.2. The rates p, q and density g G (w™™, uj"^"-^) have the fol- 
lowing properties. 

• The jump rate functions p and q satisfy assumptions (|2.6p . (|2.7p . (|2.8p . 
(j2.9p and (|2.1ip discussed in Sections \2.1\ and\ 



• n"{g) < 0. 

• Let u/) be a pair of processes in basic coupling, started from distri- 
bution p.l8[) . with second class particle Q{t). Then there exist constants 
< ao, Co < oo such that 

(2.31) p{|Q(i)| >x}<Co~ 

whenever K > agt and t is large enough. 

As mentioned, our results are valid only for asymmetric processes. The 
assumption of asymmetry is implicitly contained in Ti."{g) < 0. Symmetric pro- 
cesses have 7i(p) = 0. Exponential tail bounds for \Q{t)\ that imply assumption 
(|2.3ip hold automatically if the rates p, q have bounded increments because the 
rates for Q come from these increments of p and q. Here is the main result. 

Theorem 2.3. Let Assumvtions [K7\ and \2.2\ hold for density g. Let the pro- 
cesses {uj_~{t), Lo{t)) evolve in basic coupling with initial distribution ()2.18p and 
letQ{t) be the position of the second class particle between Lo^ (t) anduAt). Then 
there is a constant Ci = C'i{g) G (0, oo) such that for all 1 < m < 3, 

1 ,. . E\Q(t) -VHl"" ,. BIQU) - Vetl"" Ci 
2.32 — < liminf '^^ ; ^— < limsup '^^ ; ^— < 

^ ' Ci t^oo i2m/3 - ^2m/3 3 - TO 

Superdiffusivity of the second class particle is best seen with the choice 
TO = 2: the variance of its position is of order t^^^. Next some corollaries. 
Notation [X\ stands for the lower integer part of X. 

Corollary 2.4 (Current variance). Under Assumvtions [K7\ and lKM there is a 
constant Ci — Ci{g) > 0, such that 

1 Var^(V^^ Var^(V^tjW) ^ 

c; < i^sj^f — — ^ 'TJ^'p — — < 

This follows from (j2.19p with the choice to = 1. 
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Corollary 2.5 (Law of Large Numbers for the second class particle). Under 
A ssumvtions WHS and \2.2l the Weak Law of Large Numbers holds in a density-g 
stationary process: 

(2.33) ^ 4 v. 

If the rates p and q have bounded increments, then almost sure convergence also 
holds in (j2.33p (Strong Law of Large Numbers). 



The Weak Law is a simple consequence of Theorem 12.31 The Strong Law 
will be proved in Section [5l 

Corollary 2.6 (Dependence of current on the initial configuration). Under 
Assumptions [2J] and \2.Sl for any V R and a > 1/3 the following limit holds 
in the L^ sense for a density-g stationary process: 

(2 34) lim ^LVtJ W ~ ^[Vtj-Lvgtj (0) - ^(^(g) " g^'(g)) ^ ^ 



Recall that 

f 

i=[Vt\-lVet]+l 

0, if V = 

lVtl-lVH\ 
i=l 

only depends on a finite segment of the initial configuration. Limit (|2.34p shows 
that on the diffusive time scale i^/^ only fluctuations from the initial distribution 
are visible: these fluctuations are translated rigidly at the characteristic speed 
V^. Proof of (|2.34p follows by translating ftn/tj (i) — ^[ytj-[\/etj (0) to h^yet\ {t)~ 
ho{0) = hivi!t\{t) and by applying Corollary El From ^TMi . ([2351) and the 
i.i.d. initial {uji} follow a limit for the variance and a central limit theorem 
(CLT), which we record in our final corollary. 

Corollary 2.7 (Central Limit Theorem for the current). Under Assumptions 
\2.1\ and \2.2[ for any V €R in a density-g stationary process 



(2.35) hyvt\^yv^t\{Q)^ < 



(2.36) lim ^"'•'^^LV^J^^^^ ^ Var^(c.) ■\V^-V\=: 

t-*oo t 

and the Central Limit Theorem also holds: the centered and normalized height 
h^ytj {t)/^/t ■ D converges in distribution to a standard normal. 

For ASEP the CLT, the limiting variance (|2.36p and the appearance of initial 
fluctuations on the diffusive scale were proved by P. A. Ferrari and L. R. G. 
Pontes [18] . For convex rate zero range and bricklayers processes Corollary 12.71 
was proved by M. Balazs [S]. 
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Remark on the convex case 

Our results and proofs work in the analogous way in the case where the flux is 
convex and the corresponding microscopic convexity is assumed. 

2.8 Two examples that satisfy microscopic concavity 

Presently we have verified all the hypotheses of Theorem 12.31 for two classes of 
processes. 

2.8.1 The asymmetric simple exclusion process 

The asymmetric simple exclusion process (ASEP) was the first example de- 
scribed in Section 12.21 It has two parameters < p ^ q < 1 such that p + q = 1 . 
To be specific let us take p > q so that on average particles prefer to drift to the 
right. The invariant measure jj,^ is the Bernoulli distribution with parameter 
< p < 1, while /x^ is concentrated on zero for any g. The hydrodynamic flux 
is strictly concave: 7i(p) = {p — 9)^(1 — g)- 

The detailed construction of the processes y{t) and z{t) needed for Assump- 
tion [211] can be found in [12]. Here it is in a nutshell. 

Given the background process {r]{-), w(-)) and the second class particles 
{Xm{-)} between them, the processes ?/(•) and z(-) are nearest-neighbor ran- 
dom walks on the labels {m} with rates p and q. Walk y{-) has bias to the 
left (rate p to the left, rate q to the right) and walk z{-) has bias to the right 
(rate p to the right, rate q to the left). Their jumps are restricted so that jumps 
between labels m and m + 1 are permitted only when Xm and Xm+i are adja- 
cent. The clocks governing these jumps are coupled so that the ordering y < z 
is preserved. 

Since a second class particle in ASEP is bounded by a rate one Poisson 
process, (|2.31|) holds. 

We gave an earlier proof of Theorem 12.31 for ASEP in [13] . The present 
general proof evolved from that earlier one. 

2.8.2 Totally asymmetric zero range process with jump rates that 
increase with exponentially decaying slope 

This class is completely new in the sequence of models for which i^/'^-scaling of 
current fluctuations have been verified. Models in this class have a richer behav- 
ior than cither ASEP or the totally asymmetric zero range process (TAZRP) 
with constant rate. As explained in Section 12. 2i in a TAZRP one particle is 
moved from site i to site i -I- 1 at rate f{uji), and no particle jumps to the 
left (our convention for total asymmetry is p = I — q — 1). The jump rate 
/ : Z+ M+ is nondecreasing, /(O) = 0, and f{z) > for z > 0. Assume 
further that / is concave. 

As we shall see later in Section[6l one aspect of microscopic concavity, namely 
the ordering of second class particles, can be achieved for any TAZRP with a 
nondecreasing concave jump rate. Indeed, up to Lemma 16.21 in Section [6] we 
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only use monotonicity and concavity of the rates /. Thus for concave TAZRP 
only the tail control (|2.28p - (|2.29p of the label processes remains to be provided. 
For this part we currently need a stronger hypothesis, detailed in the next 
assumption. 

Assumption 2.8. Let p = 1 — q = 1. Assume the jump rate function f of a 
totally asymmetric zero range process has these properties: 

• /(0) = 0</(l), 

• / is nondecreasing: f{z + 1) > f{z), 

• f is concave with an exponentially decreasing slope: there is an < r < 1 
such that for each z > 1 such that f{z) — f{z — 1) > 0, 

(,,37) ^ 

f{z)- f{z-i) 

The case where f becomes constant above some zq is included. 



Theorem 2.9. Under Assumption \2.8l a stationary totally asymmetric zero 
range process satisfies the conclusions of Theorem \2.3\. and the conclusions of 
Corollaries \EM WM and\K^ 

A class of examples of rates that satisfy Assumption 12.81 are 

f{z) = 1 - exp(-/3z''), /3 > 0, I? > 1. 

Another example is the most basic, constant rate TAZRP with f{z) = l{z > 0}. 
For this last case a proof has already been given in [7] . 

To prove Theorem 12.91 we need to check Assumptions 12.11 and 12.21 of Theo- 
rem [531 The construction of the label processes y{t) and z{t) and verification 
of Assumption 12.11 are done in Section [Sj Assumption 12.21 requires only a few 
comments. The properties of the rates required in the first bullet of Assump- 
tion 12.21 are straightforward. Since / is concave and cannot be linear due to 
(|2.37p . Proposition lC.il in the appendix implies that 'H"{g) < for each g > 0. 
Concavity of / implies bounded jump rates for the second class particle Q(t), 
hence a simple Poisson bound gives (|2.3ip . 

The remainder of the paper is devoted to proofs. The next two sections 
prove Theorem 12. 3[ after that we prove the Strong Law for the second class 
particle, and then we return to finish the proof of Theorem l2.9l 



3 Upper bound of the main theorem 

In this section we prove the upper bound of (|2.32p . Density g is fixed. Let 
A G (p, p — 7o) and apply Assumption 12. ll with constant sequences gi = g and 
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Xi = X for all i G li. Notations P, E, Var, Cov will refer to the coupled four- 
process evolution described in Assumption l2.1i while P^, E^, Var^, Cov^ will 
refer to a density g stationary process. Abbreviate 

(3.1) ^{t):^E\Q{t)~[VH\\. 

The requirement that {u/~ , lu) obey the basic coupling was included in Assump- 
tion [TT] Consequently ^'(i) is the m — 1 expectation of (|2.32p . 

The following lemma does the main work towards the upper bound. 

Lemma 3.1. There exist positive constants ai, a2, to such that for each t > to 
and integer u such that ai\Jt < u < ait, 

(3.2) P{Q{t) > lVSt\ +u}< Cs {^{t) + u}+ Ci^. 



Proof. We start with an integer m > 0, and write 

(3.3) P{Q(i) > [VH\ +u}< P{y{t) >k} + P{Xkit) > Qit) > [V^t\ + u}. 

The event {Xk{t) > lV^t\ + u} implies that among the X^'s at most particles 
Xi, . . . , Xk^i have passed the path (s([F^iJ + u) + 1/2)^^^^^^ from right to 
left. Each such passing decreases ^"yetj+„(i) — h^ygt^^u^t) by one (recall the 
statement around (12.101) 1. Hence we can bound the probability in (|3.3[) by 

P{y{t) >k} + P{/if^,,j+„(t) - hly.^^^Jt) > -k}. 

We introduce two more processes: r]°^ is a stationary process started with 
initial data r]^^{0) — 77i(0) for i ^ 0, while 7?o''(0) is fi^ distributed independently 
of everything, uf"^ is a stationary process started with UJ°'^{Q) = uji{Q) for i ^ 0, 
and Wg''(0) is fi^ distributed independently of everything. Include these in the 
basic coupling of (ry, ui) and write 

+ f^fvetj+uW ^ ^|VetJ+«(^) 

Basic coupling implies 

htv^tl+uit) - hl^vltl+uit) < l^o(O) -^S'(0)l < ko(0)| + kS'(0)l 
and - /.j'^.^j+JO < IC(0) - '?o(0)| < |C(0)l + l^o(0)|. 

We bound the stationary expectations using (|2.16p . (|2.15p and Taylor's formula: 

= n{g)t - {[VH\ + u)g - H{X)t + ([T^^iJ + u)X 

< t{H{g)-H{X)+H'{g){X-g))+u{X-g)+Ci 

< -ln"{g)(g - A)2 + u(X - g) + C2t{g - Xf + Ci. 
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H. can be differentiated arbitrarily many times, as we show in Section [C] of 
the Appendix. Constant Ci above bounds errors from discarded integer parts. 
Recall that tilde stands for the centered random variable. Collecting terms we 
continue from (13.311 as follows. 



P{Q(t) > IVH\ + u} 

< P{y{t) > k} 

+ P{hfy,,^^.^it) - hf;,,^^Jt) >-k+ \n"{g){Q ~ \f + u{q - A) 
- C^tiQ \f ~C,- |ryo(0)| - |C(0)| - |c^o(0)| - \iol\m 

< P{y{t) > k} 

+ nht;iti+uit)-hlv\ti+uit) > \n"{e){e-x? + ^(^^- A)} 

+ P{ho(0)| + |C(0)| + K(0)| + |a;°»| 

>-k+ ^{q - A) - C2t{g - A)3 - Ci}. 

From now on we use the specific assumption Ti." (g) < 0. We maximize the terms 
on the right-hand side of the probability of ft,'s by the choice 

—u 



g — \ = 



2tn"{g)' 

To stay within the range of densities covered by Assumption 12. II we must ensure 
that X > g — jQ. So we introduce a small constant ai > and restrict our 
calculations to the case u < ait. Then 

P{Q(i) > IVH\ + u} < P{y{t) > k} 

+ P{|r7o(0)| + |C(0)l + K(0)| + 1^.^0)1 

^_ , g2 ±_rx 

An"[g) ■ t ^ W'igr ■ 

Now we set 

^8H"{g) 

and assume a2Vi < u < ait for a possibly smaller ai and a large enough a2- 
That allows us to unify the right-hand side of the inequality in the last line. 
Thus for all large u and t with a2\ft <u < a\t 

P{Q(t) > \VH\+u} < P{y{t) > • yj} 

+ P{\mm + IC(o)l + ko(o)| + IcS'^Wl > cA. 
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Assumption 12.281 allows us to bound the first probability on the right by 
CiP /u^ (take ^ — g — X). Apply Chebyshev's inequality on the second line and 
Markov's inequality on the third one: 

P{Q(i) > [VH\+u} 

< 128^^^^{Var^(/ifv.,j+J0) + Var^(/.j';\,j^JO)} 

+ €4^. 

The term C^t/u^ was subsumed under C/^fi /u^ due to the condition u < ait. 
The variances here are taken under the stationary distributions of the processes 
rf^ and 0;°*^. That allows us to apply p.l9p . whose right-hand side takes us 
back to the four-process coupling under measure P. Recall p.ip . 

P{Q(0 > \ VH\ + u} 

< cJ-^^-^{nQ{t) ~ [VH\ - u\+^\Q\t) - \ VH\ - u\\ + Ci^^ 

< C,-^^{E\Q{t) - [V^t\ I + E\Q^'{t) - [VH\ \ + 2u\ + C4- 

The variable Q^{t) above is the location of a single discrepancy between the 
process 77 and one started initially with ?y''"(0) = ?7(0) -l-^o- 

It remains to relate E|Q''(t) - [T^^^ijjto «'(i)7 This is where part (lOO]) of 
Assumption 12.11 is a key point. Compute now in the four- process coupling of 
Vi V~^j J ^ described in Assumption 12. II Use (|2.30p and Taylor expansion of 
Ti again: 

E|Q" W - IVH\ I < E(Q"(t) - Q{t)) + 
(3.4) =in'iX)-n\g))t + '^{t) 

< n"{g) ■ (A - g)t + Ceig - Xft + *(t) 

7/ 7/ 

= ^ + Ce- + *(t) <{\ + Ceai)u + *(t). 

The last inequality used u < ait. Substitute this back into the previous dis- 
play and rename constants. This finishes the proof of p.2p and completes the 
Lemma. □ 

Completely analogous arguments lead to the same upper bound for the lower 
tail of Q{t), and together we get the following bound on the tail of the absolute 
deviation, still for a2\/t < u < ait: 

P{\Qit)~lV^t\\>u}<c /'^']^^\ ^it)+u}+C4^. 
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Next we relax the restriction to integral u and the upper limit on it: 

Lemma 3.2. There are positive constants a2, such that for all t > to and 
all real u > OL2\ft, 

V{\Q{t) - \VH\\ >u}< c f^'^f'\ ^{t)^u)+cS^. 

Proof. Any u > 1 is less than twice its integer part. Hence by simply increasing 
the constants Ci, for all large t and all real u G {a2^/t^ ctii), 

(3.5) V{\Q{t)^\VH\\>u}<c f'^"[^^\ ^{t)+u}+C^^^. 

Recall ([23T|) . When a i < ao + 2\Ve\ + 2, assume aii < u < {ao + 2\Ve\ + 2)t. 
Then Q;2-\/i < u ■ ai/{ao + 2|y^| + 2) < ait for large enough t, and (|3.5p still 
holds for u replaced by u ■ ai/(ao + 2|V"^| + 2): 



^{\Q{t)-\yH\\>u}<v[\Q{t)~\VH\\ 



ai 

> u 



ao + 2ye +2, 



<C,*^^^^^{^it)+u}+cJ^ 

via modifying the constants by factors of ai/(Q:o + 2\V^\ + 2). 

Finally, when u> (ao + 2|y«^| + 2)t, the fact that u - \ \VH\ \ > aot allows 
us to use (|2.3ip : 

P{\Q(t) IVH\\ >u}< P{\Q{t)\ >u-\[VH\\} 

- \u-\\Vet\\Y - 

Combining the above cases we get the statement for all u > ol2^- CH 

Proof of the upper bound of Theorem \2.3[ We now fixr>0, l<m<3, and 
write 



^{\Qit) - [VH\r) 

oo 

= J P{\Q{t)-lV<^t\r>v}dv 



4 — TO 3 — TO 

First choose m = 1 and r large enough to get ^'(i) < Ct^/^. Then insert 
this bound back into the last line of the display to get the bound for general 
1 < TO < 3. □ 
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4 Lower bound of the main theorem 



In this section we prove the lower bound of (|2.32p . Density g is fixed again, 
and A € (f? — 7o, p) is a varying auxihary density. We let the jointly defined 
four processes (77, ry"*", w^, lo) be exactly as defined in the upper bound proof of 
Section [3l namely, as given by Assumption 12.11 with constant densities \i = A 
and Qi = g. The initial distribution of (77, u) is /i-- of (|2.24p . Two second class 
particles start from the origin: between processes 77 and 77"*", and Q between 
processes u/~ and lu. The quantity of primary interest is abbreviated, as before, 
by*(t)=E|Q(t)-Ll/^tJ|. 

To prove the lower bound of (|2.32p it suffices, by Jensen's inequality, to 
prove the case m ^ 1. This means showing that ^'(i) > Ct^/"^ for large t and a 
constant C > 0. 

4.1 Perturbing a segment initially 

For this proof we need to introduce another coupled system and invoke Assump- 
tion [^T] once more. By concavity of the flux characteristic speeds — H^g) 
and — TL'{\) satisfy < V^. Throughout this section u > denotes a 
fixed positive integer, and 

n = lVh\ - IVH\ + u. 

Recall definitions (|2.22p and (j2.23p of the single-site coupling measures. Let 
(^(■): C(')) be a pair of processes that obeys the basic coupling, and whose initial 
distribution is the product measure 

i<—n i— — n —n<i<Q Q<i 

This initial measure complies with the pattern in (|2.24p . but translated 77 sites 
to the left so that the distribution at site ~n instead of the origin. A 

few points about this initial state: f(0) has the stationary density-A product 
distribution except at site —77, where it is ^^-distributed. C(0) has the product 
distribution with marginals /i^, except at sites {— n -I- 1, . . . , 0} where the pa- 
rameter g switches to A, and at site —n where it has distribution + 1. At 
sites -77 < 7 < n^-^ forces ^.^(0) = C»(0). 

We add a second class particle to the process ^(•), start it at site —77 and 
denote its position at time t by Let ^^(t) := ^{t) + SQ-,^(^ty 

As described in Section 12.61 the C ~ C second class particles are labeled and 
their ordered positions denoted by {Xm{t)}. The labeling is chosen to satisfy 
initially 

(4.1) • • • < a:_i(o) < a:o(o) = -77 < o < a:i(o) < X2{o) < . . . 

Thus initially Xo{0) — —n = Q(^")(0). We invoke Assumption 12.11 to have a 
label process z{t) with tail bound ((2^ such that Q'^^'^^t) = X^(t){t). (Here ^ 
plays the role of rj and C, plays the role of cj of Assumption 12. II) . 
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As before, the heights (or currents, recaU (|2.5p ) of the processes ^(•) and 
(■(•) are denoted by h^^y^j and ^^^^j, respectively. The first observation is that 
gives one-sided control over the difference of these currents. 

Lemma 4.1. For any V E M. 

Q^-'^^t) < lVt\ implies hfy^^ {t) - h^y^^ {t) < -z{t). 

Proof. Recall again, from (|2.5p and the statement around (|2.10p . that the height 
difference ^'^y^j (t) — ^^y^j (t) equals the net number of second class particle 
passings of the path (s[FiJ + 1/2)^^^^^^ from left to right. That is, each left- 
to-right passing increases ^[y^j (t) — h^'^y^^ (t) while each right-to-left passing 
decreases it. 

Suppose z{t) < 0. Then gUJ and X^(f)(t) = < L^^J imply that 

only those second class particles with labels z{t) + 1, z{t) + 2, . . . , could have 
passed the path (s[V^tJ + 1/2)^^^^^^^ from left to right. The claim follows. 

If z{t) > 0, then = Q^""^*) < L^^J implies that at least those 

second class particles with labels 1, 2, ... , z{t) have crossed the path (s[V^iJ + 
l/2)p^^^j^ from right to left. Again the claim follows. □ 

Let S(-) be a process started from the product distribution ( /i^) (g) 

1). The next lemma compares the distributions of ^ and Q_. No coupling of C 
and Q is proposed or required. 

Lemma 4.2. There exist constants 7 = 7(g) > and Ci{g) < 00 such that for 
all X Cz (f? — 7, g) and all events A the following inequality holds: 

P{Ce ^} < P{^e A}5 •exp{Ci(£-)n(f?- A)2}. 

Proof. We use the Cauchy-Schwarz inequality below to perform a change of 
measure on the distribution of the C process. First we condition on the initial 
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^-configuration at sites {—n + 1, —n + 2, . . . , — 1, 0}. 

PKeA}^ P{Ce A|C-„+i(0) = z-„+i,...,Co(0) = zo} 





[ n /^'(^'). 



n 



^= n+1 



^[ E [P{CeA|C-„+i(0) = z_„+i,...,Co(0) = zo}]' n ^'(^^) 



•S-n.+ l,---)^0 



-i— — n+1 



[ E n 



1 



< 



z_„4.i,...,zo i=-n+l 

[ E P{Ce A|C-n+i(0) = z-„+i, 



i— — n+1 



z_„+i,...,zo 



z_Ti + i,...,zo i——n+l ^ ' 

Tlie last inequality came from dropping the square. For the last equality 
note that the distributions of the initial configurations {0)^(0)} and {Ci(0)} are 
product-form and agree outside the interval — n + 2,...,— 1, 0}. Thus 

conditioned on the initial values in — n + 2,...,— 1, 0} these processes 

have identical conditional probabilities. 

To complete the proof we bound the last factor in brackets. Recall formulas 
(|2.12p and (|2.13p for the state sum and the site-marginals. Without the power 
1/2 the factor in brackets equals 



E 

n + l,---ZO 



n 

— n+l 



= - 



In the appendix we show that logZ(6') and B(^q) are infinitely differentiable. Let 
e — 9{q) — 9{\). By local Lipschitz continuity of the function 9{q), the interval 
{9{\) — e, 9{\) + e) is in {9_, 9) with a small enough choice of 7. There exists 
some 9 G {9{\) - e, d{X) + e) such that 



log( 



Z{29iX)-9ig))Z{9ig))\ 



Z{9{\W 



= log Z{9{\) - e) + log Z{9{\) + e) 
- 2\og Z{9{\)) 

= \^\ogZ{9)e^<C,{e)-{e~^? 
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Thus we get the bound 

fZ{26{\)-e{Q))Z{6{Q))\^ , ^ , ^,21 

( Z{d(\)Y ) <exp{Cl(^.)•n(^.-A)2}. □ 

4.2 Completion of the proof of the lower bound 

The gist of the proof is to get upper bounds on the complementary probabihties 
P{(5("")(t) > \yHW and P{Q(-")(t) < \VHW. As stated u is an arbitrary 
but fixed positive integer and n — [V^'^iJ — L^^^J + ^• 

Lemma 4.3. 

P{Q(-")(t)>Ll-^tJ}<M + ^Hfc^+2_ 



Proof. Distribution-wise the system is a translate of (ry, 77+, Q''), 

and so 

P{Q(-") (t) > IVH\ } = P{g(-") {t) + n- lVh\ > u} 

E{\Q'Kt)-[v't\\) 



PIQ^W - iv^t\ >u}< 



u 



^ E(|Q"(t)-Qft)|) ^ E{\Q{t)~lVn\\) ^ IVH\~IVH\ ^ 
~ u u u 

Use (I2.30p precisely as was done in (|3.4I) to conclude that the first term equals 

u-iE(Q"(t) - Q{t)) = u-H{H'{\) - H'{q)) = -u-^H"{v)t{g - A) 

for some v E (A, g). The second term is 'i>{t)/u, and the third term is similarly 
estimated by —u~^H" (y)t{Q — A) + 2/it, the last part coming from discarded 
integer parts. Setting C2 := 2 max —H"{v) finishes the proof. □ 

Notice that H"{g) < was crucial in the previous proof, as well as in the 
following lemma, and the final proof thereafter. These points show where the 
proof fails for symmetric systems ~ recall that these would have lower-order 
current fluctuations on the characteristics. 

Lemma 4.4. Let K be any number such that < K < —^tH"{g){g — A)^. 
Then for small enough 7 > 0, large enough t, and A G {g — J, g), 



^tH"{g){g- \)^ -K 



3 

C4 



\tH"{g)[g-\Y-C^t{g-\f-g 



■ e 



6 



Var'^(7?o)*W , C^t{g-\) , Q 



5 



if2/4 K"^ A'-4|A| 

+ Ct'^-'^-f^''-^K-''. 
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Proof. Lemma 23] leads to 

(4.2) < P{-z{t) > K/A] 

(4.3) + P{/if^,,j {t)<K + t{H{\) - XH'ig))} 

(4.4) + P{hl^,,^ (t) > 3K/4 + t{HiX) - XH'ig))}. 

To boimd (|4.2p we use the assumed distribution bound (|2.28p on z{t) and get 

P{-z{t) > K/A} < Ct^-^-f^''-^K-''. 
Apply Lemma 14.21 to line (|4.3p to bound it by the probability of the process S: 

(ESI) < [P{/ify.tj(t) <K + t{H(X) - XH\g))}]^ ■e^'''^s-^'>\ 

As in the proof of Lemma 13.11 we switch to stationary processes to get precise 
bounds: 

+ [mfy,,^ it) - t{H{g) - gH'{g))] + t{H{g) - gH'{g)) 
> hf,,^ (t) - P_„(0)| - Icvlim - \g\ + t{H{g) gH\g)). 

After the equality sign, the absolute value of the first term in brackets is not 
larger than |tD_„(0)—(jjl'^„ (0)1 < |tD_„(0)| + |ci;!l'^„(0)|. The second term in brackets 
is between —\g\ and \g\ due to the integer part in [V^^tJ. Consequently 

K(y,^^{t)<K + t{H{X)^XH\g)) 

implies 

h-yy,,^ (t) - p_„(0)| - livlim <K + t[HiX) H{g) + H'{g){g A)] + |^,| 

<K+ \tH"{g){g - X)^ + C^t{g - A)^ + \g\. 

Then, we cut the event into two parts according to the value of |a)_,i(0)| + 
l'^-'ri(0)l ^''^d we use p.l9p to bound the variance of Var[/i'iygjj {t)] by the func- 
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tion 

+ [P{P_„(0)| + Ic.-jO)! > -\tH"{Q){e Xf ~ C^tig - Xf - \g\}] ' 

- A)2 - if 



e 



-^tH"{g){g-Xr~C3t{g~Xr~\g\ 

Now we turn to (|4.4|) . To reduce ft-^ye^j to the current of the density-A equi- 

hbrium process ^^y^^j and to get rid of the integer part errors we argue as 
before. 

"-[Vet] — "-[Vet] + Velvet] ""[Vetji 

+ [E^/i;'y,,j - t{H{X) - XH'ig))] + t{H{X) - XH'{g)). 

hfy^^^it) differs by at most |e-„(0) - r/lim < |C_„(0)| + \v-im from 
/ij'ygjj (i)). Taking integer parts again into account, giving another error term 
|A|, Hne (|4.4p is bounded from above by 

^{hlvnlit) + ie-n(0)l + \rtl{0)\ + \X\ > SK/a]. 

Then, we cut the event into two parts and use Markov's inequahty in the second 
one: 



O) < p 



{h'ff.tiit) > ^/2} +P{|e-„(0)| + Ivlim > K/4- 



< 



X~4|A|' 

We can use p. 191) again to continue with 

Var\^o)n\Q''it)-[VH\\) , C5 



< 



K^/A K-A\X\' 

Repeating the first two steps of calculation p.4p we can write 

E{\Q^{t) - IVH\ I) < E{\Q^{t) - Q(i)l) + E(|g(t) - [VH\ |) 
<Ct{g- XJ + 'i'it). 
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So, we finally get 

K-A\X\' 

Proof of the lower bound of Theorem \2.3[ As observed in the beginning of this 
Section, it suffices to prove that 

(4.5) liminft"2/3^(^^^ Q_ 

t — >oo 

In the last two lemmas take 

u = [M^/^l , g-X = bt-^'^, and K = bt^'^ , 

where h and b are large, in particular b large enough to have b < ~^H"{g)b^ so 
that K satisfies the assumption of Lemma 14.41 Then 

n = lVh\ ~ IVH\ + u 

< {H' {X) - H' {g))t + u + 2 

= -H"{q){q -X)t + u + Cjtig -Xf+2 

< {-H"{g)b + h)t^/^ + Cyfe^i^ + 3 

for large enough t. With these definitions we can simplify the outcomes of 
Lemma 14.31 and Lemma 14.41 to the inequalities 

(4.6) P{Q(-") it) > iV^tl } < + ^ + ^ 
and 

(4.7) P{Q(-")(,) < L^.^J} < '^\c^ + ^ + ^ + Cb^-^ 

The new constant C depends on b and h. 

The lower bound (|4.5p now follows because the left-hand sides of (|4.6p - (|4.7p 
add up to 1 for each fixed t, while we can fix b large enough and then h large 
enough so that C26//1 + Ce/b + Cb"-^ < 1 (recall k < 3). Then t-^/^^!{t) must 
have a positive lower bound for all large enough t. This completes the proof of 
Theorem El □ 



5 Strong Law of Large Numbers for the second 
class particle 

This section proves the Strong Law of Large Numbers (CoroUarv 12. 5p . We 
assume that the jump rates of the second class particle are bounded, i.e., 
(5.1) 

p{y+l,z)-p{y,z), p{y,z)~p{y,z + l)]^ V..-'" < y, . < c.— . 
<liy, z + l) - q{y, z), q{y, z) - q{y + 1, z)\ 
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This means that the second class particle has at most rate C to jump to the 
right and to the left, respectively, implying that starting at any time i, it can 
be bounded by rate C Poisson processes that start from its position Q{t). 



Proof of Corollary \2.5[ Let e, S > 0. Define the events 



An ■ — 



,1+5 



>e/2} 



for n £ N. Then, Markov's inequality and Theorem 12.31 imply, for 1 < to < 3 
and large n, 

P{A„} = P{\Q{n^+^) ~ F^ni+'^l" > (e/2)"n(i+'')'"} 
1 



< 



< 



(£/2)™n(i+'5)™ 
Ci 



F,[\Q{n^+^)-VSn^+^r] 
1 



(3-m)(e/2)™ „rn(l+'5)/3 ' 

which is summable if (1 + d)m > 3. Here d can be chosen arbitrarily small by 
taking m close to 3. By the Borel-Cantelli Lemma there exists a.s. no G N such 
that 



(5.2) 



Vn > no 



,l+<5 



ye 



< e/2. 



Using this we show that a.s. there exists ni e N such that 



(5.3) 



Q{t) 



< £ for all real t > 



Let n > Uf) and suppose there exists some t e [n^^^ , (n + 1)^^*^) such that 
fails: \Q{t) - V<^t\ > et. Together with we have, if n is large, 

(5.4) 

\Q(t) - Q{n^+')\ > \Q{t) ~ VH\ - \Q{n'+') - V^n'+'\ - \V^t - V'n'+'\ 



>et~e/2-n 



1+5 



\v'\{t 



,1+5 



> —n 
~ 4 



1+5 



The jump rates (|5.ip (both left and right) of Q are bounded by C. However, 
the event (|5.4p implies that at least [|n^+''J many left jumps or this many right 
jumps happen in the time interval [n^^^ , (n+ 1)^+*). For large n, the length of 
this interval is smaller than 2(1 + S)n^ . Let N{-) be a rate C Poisson process. 
Then for large n the probability of the event (|5.4p is bounded from above by 

2P{N{2{1 + 6)n') > |nl+*} < 2P{e^(2(l+5)r/) > ^e/4-n^ + '^ 



< 2e-^/4-"'+' 
= 2e-^/4-"'+' 



gjgJV(2(l+5)«'')j 
. (c-l)2C(l+5)n'' 



This quantity is summable over n, so the Borel Cantelli Lemma implies that 
a.s. ()5.3p holds eventually. Since this is true for each e > 0, the Strong Law of 
Large Numbers holds. □ 
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6 Microscopic concavity for a class of totally 
asymmetric concave exponential zero range 
processes 

In this section we verify that Assumption 1 2 . 1 1 can be satisfied under Assumption 
12.81 and thereby complete the proof of Theorem 12.91 

The task is to construct the processes y{t) and z{t) with the requisite 
properties. First let the processes (??(•), evolve in the basic coupling so 

that rii{t) < uJi{t) for all i S Z and t > 0. We consider as a background 
process this pair with the labeled and ordered u — rj second class particles 

• • • < X^2{t) < < Xo{t) < Xl{t) < X2{t) < ■ ■ : 

At each time t > this background induces a partition {Aii{t)} of the label 
space Z into intervals indexed by sites i (li 1^, with partition intervals given by 

M^it) := {m : X,n{t)=i}. 

(For simplicity we assumed infinitely many second class particles in both direc- 
tions, but no problem arises in case we only have finitely many of them.) Mi{t) 
contains the labels of the second class particles that reside at site i at time t, and 
can be empty. The labels of the second class particles that are at the same site 
as the one labeled m form the set Mx„,{t){t) = : {a"(i), a'^it) + 1, • • • , b'^it)}. 
The processes a™{t) and are always well-defined and satisfy a^"'{t) <m< 

lf^{t). 

Let us clarify these notions by discussing the ways in which a"^{t) and 6™(i) 
can change. 

• A second class particle jumps from site Xm{t—) — 1 to site X„i{t—). Then 
this one necessarily has label a"^{t—) — l, and it becomes the lowest labeled 
one at site Xm{t—) = Xm{t) after the jump. Hence a™{t) — a™{t—) — 1. 

• A second class particle, different from X^, jumps from site Xm{t—) to 
site Xm{t—) -1-1. Then this one is necessarily labeled b"^{t—), and it leaves 
the site X„,{t-), hence = b"'{t-) - 1. 

• The second class particle Xm is the highest labeled on its site, that is, m = 
b'^{t—), and it jumps to site Xmit—) + 1. Then this particle becomes the 
lowest labeled in the set Mxm(t-)+i — ■^x„(t)j hence a™(t) = m. In this 
case can be computed from 6™ (t) — a™ (t)-}-l — iLjx„,{t){t)-'ilx^{t){t), 
the number of second class particles at the site of X^, after the jump. 

We fix initially y(0) = z{0) = 0. The evolution of {y, z) is superimposed 
on the background evolution (77, w, {Xm}) following the general rule below: 
Immediately after every move of the background process that involves the site 
where y resides before this move, y picks a new value from the labels on the 
site where it resides after the move. Thus y itself jumps only within partition 
intervals Mi. But y joins a new partition interval whenever it is the highest 
X-label on its site and its "carrier" particle Xy is forced to move to the next 
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site on the right. This is the situation when y{t~) — and at time t 

an oj — rj move from this site happens. (Recall that the choice of X-particle to 
move is determined by rule (|2.25p . In the present case there is only one type of 
Lo — rj move: the highest label from a site moves to the next site on the right.) 
All this works for z in exactly the same way. 

Next we specify the probabilities that y and z use to refresh their values. 
When y and z reside at separate sites, they refresh independently. When they 
are together in the same partition interval, they use the joint distribution in the 
third bullet below. 

• Whenever any change occurs in either or 77 at site ) and, as a 

result of the jump, a^*^*~)(t) 7^ a^'*^'(t), that is, y{t—) and z{t~) belong 
to different parts after the jump then, independently of everything else, 
(6.1) 

a^(*-)(t) -.-ithpr 

/(^x„(t_)(t)(0) - /fe„(t-)(t)(i)) 
/(^x„(,_,(t)(0) - /(^x„(,_,(t)(0) 



v{t) 



when the denominator is non-zero, and yit) : — a^'* ^ [t) when the denom- 
inator is zero. 

• Whenever any change occurs in either w or 77 at site Xzi^t_-^{t~) and, as a 
result of the jump, a^(*")(t) ^ a^(*~^(i), that is, y{t~) and z{t-) belong 
to different parts after the jump then, independently of everything else, 
(6.2) 

{t) - 1 with pr ^('^^-(*-)(*)(^»-^(^^-(*-)(*)(^) + ^) 

b<'~Ht) withpr + 

/(wx,(,_,(t)W) - /('7x,(,_,(t)(0) 



z{t) 



when the denominator is non-zero, and z{t) : — b^^*^^\t) when the denom- 
inator is zero. When a;x,(t_)(t)(i) — ''lx^^t-){t){t) + li h^'^^~\t) — 1 is not 
an admissible value but in this case the probability in the first line is zero. 

Whenever any change occurs in either u; or 77 at sites Xy(^i__-^[t~) or 
X^^^^_){t—) and, as a result of the jump, a^(*~^(t) = a^(*~)(t), that is, y{t—) 
and z{t—) belong to the same part after the jump, that is, Xy(^i_){t) = 
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yit) 

z{t) 



Xz(t-){i) then, independently of everything else 
(6.3) 

^5y(t-)(t) _ 1^ 

with pr. ^^^^f""-'^*^^^^,t'^'/^'""''"^;\^!^^ 

_ /(^X„(,_)(t) W) " /(^X„(t-)(t) W - 1) 

V(*-)(t)\ 

when the denominator is non-zero, and 

(y(t),zW):=(a'^(*-)(t), 

when the denominator is zero. When wx^(t„)(t)(0 = Vx^(t-)(t)(^) + 1> 
b^^*^\t) — 1 is not an admissible value but in this case the probability in 
the first line is zero. 

The fact that the numbers on the right hand-sides are probabilities follows from 
u)i{t) > ri.i{t) on the sites i in question, and from the monotonicity and concavity 
of /. The above moves for y and z always occur within labels at a given site. 
This determines whether the particle Q{t) := Xy(^-^{t) or Q^{t) : = Xz(^t){t) is 
the one to jump if the next move out of the site is an a; — 77 move. 

We prove that the above construction has the properties required in As- 
sumption 12.11 



Lemma 6.1. The pair (w , ttl) : = ^i^) obeys basic coupling, as does 

the pair (77, 77+) : = (77, 77 -|- S^J- 

Proof. We write the proof for (w^, ui). We need to show that, given the config- 
uration (j], u, {Xm}, y), the jump rates of uj) are the ones prescribed in 
basic coupHng (Section [2.3p and by \2.2\ . Leftward jumps of type (|2.3|) do not 
happen in the system under discussion. Since the jump rate function p depends 
only on its first argument, jumps out of sites i ^ Q happen for cj^ and a; with 
the same rate p{uj~ , ^i+i) = /(w~) = f{i^i) — p{^i, ^i+i)- The only point to 
consider is jumps out of site i = Q. 
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Since the last time any change occurred at site i, y chose values according to 
(|6.ip or (|6.3p . Notice that l|6.1[) and (|6.3[) give the same marginal probabilities 
for this choice. Hence 

(6.4) y took on value with probability — ^ — f^J^^ 

/K) - fivi) 

and 

/(w,) - fiiu, ~ 1) 



(6.5) y took on value 6^ with probability 



as given in (|6.ip . or y took on value in the case f{uji) = fiVi)- According 
to the basic coupling of rj and lo, the following jumps can occur over the edge 

• With rate p(wi, aji+i)-p(?7i, ry^+i) = f (uji) - f (rii) , when positive, w jumps 
without rj. The highest labeled second class particle, Xtv jumps from site 
i to site i + l. 



— With probability (|6.5p Xj, = Q jumps with X^v ■ In this case 
since the difference Q disappears from site i. Also, 

^i+iit-) = ^t+iit~) = Wi+i(t) - 1 = Wi:^_i(t), 

since the difference Q appears at site i+l. So in this case lu undergoes 
a jump but lo~ does not, and the rate is 



— With probability (|6.4p = Q does not jump with Xf,y, since it has 
label and not (this probability is zero if coi = rji + 1). In this 
case and lj perform the same jump and it occurs with rate 

-/(,.)! ■ ^^;,^;,^y =/(.r) -/(.). 

• With rate p(?7i, ?7i+i) = fiVi), both and jump over the edge (?, z + 1). 
No change occurs in the lu — rj particles, hence no change occurs in Q. 
This implies that the process jumps as well. 

Summarizing we see that the rate for (w~, ui) to jump together over (z, z + 1) is 
f{ui^), and the rate for a; to jump without ui" is /(wi) — f{uJ^)- This is exactly 
what basic coupling requires. 

A very similar argument can be repeated for (77, 77+). □ 
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Lemma 6.2. Inequality (j2.26p y <^ z holds in the above construction. 

Proof. Since no jump of ?/ or z moves one of them into a new partition interval, 
the only situation that can jeopardize (I2.26P is the simultaneous refreshing of y 
and z in a common partition interval. But this case is governed by step (j6.3p 
which by definition ensures that y < z. □ 

So far in this section everything is valid for a general zero range process with 
nondecreasing concave jump rate. Now we use the special convexity requirement 
(|2?37)) . With r e (0, 1) from ((2?37)) . define the geometric distribution 



(6.6) z/(to) 




TO > 
TO < 0. 



d d 

Lemma 6.3. Conditioned on the process (77, u/), the bounds y{t) < v and z{t) > 
—V hold for all t> 0. 

The proof of this lemma is achieved in three steps. 

Lemma 6.4. Let Y be a random variable with distribution v , and fix integers 
a < b and rj < lo so that uj — rj = b — a + 1. Apply the following operation to Y : 

(i) if a < Y < b, apply the probabilities from (j6.1|) [equivalently, (j6.4p and 
(|6.5p ) with parameters a, b, rj, lu to pick a new value for Y ; 

(ii) ifY<a orY > b then do not change Y . 

Then the resulting distribution v* is stochastically dominated by v . 

Proof. There is nothing to prove when b ~ a, hence we assume & > a or, 
equivalently, w — 77 = 6 — a + l>2. It is also clear that i>* (m) — v{ni) for m < a 
or TO > 6. We need to prove, in view of the distribution functions, 

m m b b 

J2 ^* W ^ H "'^^^ equivalently, ^ < ^ iy{e) 

£—a I— a £—in i—m 

for all a < TO < 6. Notice that v* gives zero weight on values a < m < 6 (if any), 
therefore the left hand-side of the second inequality equals i'*{b) for a < m < b. 
Hence the above display is proved once we show 



iy*{b) < v{b), that is, 



f{L0) - f{LO - 1) 



b 



(6-7) ^^TZTT^^TTZ^ • 2. ^^^) ^ ^^^) 

1= 



see (|6.ip . When /(w) = f{u) — 1), there is nothing to prove. Hence assume 
/(w) > fiio— 1) which by concavity implies that / has positive increments on 
{77, . . . , uj}. If & < then both sides are zero. If > then we have, by (|2.37p . 

n 

f{uj-b + £)-f{uj-b + e~i) 
' ■ /(^) - /(^ - 1) 
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for each £ < b. The first inequahty also takes into account possible v(£) = 
values for negative £'s. With this we can write 

which becomes (|6.7p via w — 77 = 6 — a+1. □ 

We repeat the lemma for z{t). 

Lemma 6.5. Let Z be a random variable of distribution —v, and fix integers 
a < b, T] < Lu so thatuj — rj = b — a+1. Operate on Z as was done for Y in Lemma 
\6.4\ but this time use the probabilities from (|6.2p with parameters a, b, r], lu. Let 
—v* be the resulting distribution. Then v* is stochastically dominated by v. 

Proof. Again, we assume 6 > a or, equivalently, uj — rj — b — a + l>2. It is also 
clear that m) = v{—m) for m < a or m > b. We need to prove 

m rn 
t=a l=a 

for all a < m < 6. Notice that —v* gives zero weight on values a < £ <b— \ (if 
any), therefore the left hand-side of the inequality equals for a < rn < 6 — 1, 
v*(\) — 1) for m = 6 — 1, and agrees to the right hand-side for m = b. Hence the 
above display is proved once we show 

> iy{-b), that is, 

see dO]). We have, by ((237)) . 

/(^)>/(z-l) 

fiv + l + b-£)-fiv + b-£) 

^ fiv + i)-fiv) 

for each £ < b. The first inequality also takes into account possible &) = 
values for positive 6's. With this we can write 

^ /('7 + 1) - /W 

which becomes (|6.8p via oj — rj = b — a+1. □ 
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Lemma 6.6. The dynamics defined by (|6.ip or (|6.2p is attractive. 



Proof. Following the same realizations of (|6.ip . we see that two copies of y{-) 
under a common environment can be coupled so that whenever they get to the 
same part A4i, they move together from that moment. The same holds for 
z(-). □ 

Proof of Lemma \6.3[ Initially y{0) = by definition, which is clearly a dis- 
tribution dominated by v of (|6.6p . Now we argue recursively: by time t the 
distribution of y{t) was a.s. only influenced by finitely many jumps of the en- 
vironment, which resulted in distributions vi, then 1^2, then v^, etc. Suppose 

d 

Uf, < v, and let v* be the distribution that would result from v by the k + 1^' 

d d d 

jump. Then Vk+i < v* hy Vk < v and Lemma 16.61 while v* < v hy Lemma 16.41 
A similar argument proves the lemma for z(-). □ 



Appendices 

A Convexity and total positivity 

This section derives a general convexity result for exponentially tilted measures. 
Let V he a, nondegenerate probability measure on R and assume that for some 
open interval X C M, 

(A.l) Y{e) = j e""^ v{dx) < 00 for all 6 el. 

For 9 £ I define the exponentially tilted measures ly^ by 

(for bounded Borcl test functions g). The nondegeneracy assumption (that v is 
not supported on a single point) and (jA.ip guarantee that 




is a finite, continuous, strictly increasing function that maps T onto a nontrivial 
open interval J. For q € J the inverse function is denoted by 9{q). 

Let ■0 be a measurable function on R, and assume (by shrinking X if neces- 
sary) that 

J\^\dv'^ <oo for all 6* el. 

Since |a;|'^ < fc!e^'^(e^^ + e^^^) for any e > and I is an open interval, it follows 
that < 00 for all fc > and 6 £j. Consequently as a function of 9 

the integral J ip dv^ has derivatives of all orders. 
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A particular case is ip{x) — x which gives the infinite difFerentiabiUty of g{9). 
Let us also note the infinite differentiability of the inverse function 6{g). Since 
g'{0) is the variance of the distribution , q'{0) > by the nondegeneracy of v, 
and so directly from the definition of the derivative 6*' (p) = 1/ g'{9{g)). Repeated 
use of basic differentiation rules produces all derivatives 0*^"^ (g). Notice that this 
argument shows a uniform lower and upper bound of g'{9), that is, Lipschitz 
continuity of both g{6) and 9{g) on bounded closed intervals. 

Define 



^ is also infinitely differentiable as a composite of two such functions. 

Theorem A.l. Assume ^ is a convex function on M. Then is convex on J . 
Assume furthermore that no linear function g{x) — aX'\-b satisfies ip = g u-a.e. 
Then {g) > for all g Cz J7 and in particular ^' is strictly convex on J . 

Proof. The proof can be reduced to the theory of total positivity. In what 
follows, citations and terminology are from Karlin's monograph |24| . The claims 
made in our Theorem I A . II follow from applying Theorem 3.5(a)-(c) from p. 285 
of [21] to the operator 

Ti:{g) = [ V'd:/''^^) = / K{g, x)^P{x) i^{dx), gej, 



where the kernel is defined by K{g, x) = Y{g{9))~^e^'^^^^ . The property of the 
kernel K that gives the result is extended total positivity (ETP) of order 3. 
This is the requirement of strict positivity on certain types of determinants of 
partials of dimensions up to 3 x 3: for all {g, x) G J' x M., 



n entries 





■ ,q\ 


= det 


\x, . 


.,xj 


l<z,j<n . 



dg^~^dx^ 



— K{g, x) 



>0 for n= 1,2, 3. 



(A.2) K 

We argue this in stages. 

We first observe that the kernel L{9, x) = Y{9)-'^e^'' on J x R is ETP of 
all orders. Recall that the Wronskian of n functions /i, . . . , /„ is the n x n 
determinant 

W[f,,...Jn]ix)^ det [ft'\^)]- 
If u is another function, the Wronskian satisfies the identity 

(A.3) W[ufi,...,ufn]{x)^u{xrW[fi,...Jn]ix). 

To justify (jA.3p . Leibniz's rule 

{uf.Y'-''=t(il'^f {l<J<n) 

k=l ^ ^ 
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implies that the matrix A = [(u/i)*--' ^•'(^)]i<i,i<n is the product of the matrices 
B=UY''-'\x)],<,,,<,, and C = [{iZlW^-'H^Mk < j}]^^^^^^^. 

By upper-triangularity detC — u(a;)". Then the corresponding determinant 
identity det(A) = det(B) • det(C) is precisely (|A.3|) . 

Now we can verify the ETP property of kernel L, utilizing (jA.SP : 



det 



dx 
= Y{9) 



TTQei^L{e, x) 



-"e"*'^l^[l, 



= det 

1 < i , J < n 
/in — li 



Ol- 



der 



= y(6l)-"e"^^ W j\ > 0. 
i=i 

To go from L{9, x) to K{g, x) = L{d{g), x), consider the 3x3 determinant 
that appears in (|A.2p , apply the chain rule and a row operation: 



K 




^XX 


Kg 


K 

QX 


K 

QXX 


Kgg 


K 

QQX 


K 

QQXX 





L 


Lx 






Lxx 








Le^g 


Lex 


Oq 






LexxOg 






Lee 


9g + LeOgg 


Leex^l 4- 


Le 


QQ ^ 


wexx&l + Lexx 


9 

Jgg 




L 


Lx 


Lxx 






L 


Lx Lxx 




Lee 


B LexOg 


LexxOg 






Le 


Lex Lexx 


> 


Lee 


9g Lega:9g 


Leexx^l 






Lee 


Leex Leexx 





The last inequality is by the ETP property of kernel L and the strict positivity 
6'g > of the derivative. The 1x1 and 2x2 determinants in (IA.2P are principal 
minors of the determinant above and are positive by the same reasoning. 

We have shown that the kernel K has the ETP property of order 3. In 
addition to ETP, Theorem 3.5 from p. 285 of [23] requires the hypotheses 



K{g, x) i'{dx) = 1 and / K{g, x)x i'{dx) — ag + b 



for some a > and 6 G M. The first one is true by virtue of the normalization 
Y{9)^^, and the second one with a — 1 and 6 = by the definition of p{9). The 
proof is now completed by an appeal to Theorem 3.5 from p. 285 of [2|]. □ 

These convexity properties can also be proved in an elementary way by 
developing suitable correlation inequalities. Such a proof is given in the note 
[TD] . We are indebted to an anonymous referee of that note for pointing out the 
connection with total positivity. 

Subsequent sections of the appendix extract from Thcorcm lA.ll consequences 
for the processes we study. 
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B Monotonicity of measures 

In this part of the appendix we show that the measures and fl^ defined in 
(|2.13p and (|2.17p . respectively, are stochastically monotone as fimctions of g. 
We start with a simple 

Lemma B.l. Fix a function if(uj) on bounded by a polynomial. Then 
E*((^(w)) is differentiable in on {9_, 9), and 

^E'Mcj)) = Cov^(^(c.), c.). 

Proof. Convergence of the series involved in E^{(p{llj)) can be verified via the 
ratio test, even after differentiating the terms. Since /i^ is the exponentially 
weighted version of for some 6*0, we have 

_d_ g d E'^"(y>(cj) ■e(^~'^«)'^) 

E^o(^H.a;.e(^-^o)") e^q . e(^-^o)^) 

Eeoe(«-eo)i^ ' ' ' [Eeoe(e-9o)c^]2 

= Cov''(^(cj), LU). 
Corollary B.2. For any 9_< 6 < 9, the state sum (|2.12p satisfies 
(B.l) ^logZW - ^ =E^.) 

(B.2) ^log^W = ^^^(e) = Var^H. 

The function g{9) is strictly increasing and maps (9, 9) onto (w™™, tj™"^). 

Proof. Everything is already covered except the last surjectivity statement. Due 
to the monotonicity and continuity one only needs to show convergence at the 
boundaries 9, 9 to w™'", w™'*'^. First let us consider the case when 9 < oo. Then 
^max _ ^ Fatou's lemma implies 

liminf Z{9) = liminf > 7 liminf ., = > = 00 

e^S B/B /(z)! e^B f{z)\ ^ /(z)! 



□ 



since for 2: > 



/(!/) 

by definition of 9 and / being nondecreasing. This shows that log Ziff) takes 
on arbitrarily large values as 6* 9. We also know that it is a smooth and 
convex function on (9_, 9) (see (|B.2P ). This implies that its derivative (IB.ip is 
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not bounded from above i.e., arbitrarily large q values can be achieved. The 
same reasoning works in case 6_> —oo for arbitrarily large negative q values. 

When 9 — oo then, regardless whether o;™'*'^ is finite or infinite, fix any 
< y < a;™'*'' and write 



(B.3) 



gie) = ^"{uj ■ l{uj > y}) + E^([cj]+ • l{uj < y}) - 'E\[uj]- ■ l{uj < y}) 
> (y + 1) • P<'{lo >y)~ E»(M- • 1{lu < y}) 

>iy + l)-iy + l)-P'ito<y)~ ^E^OH^ • ^F<^{u; < y) 
>{y + l)~{y + l)-P'{^<y)~ ^^^"m~f) ■ ^P^{u: < y) 



for a fixed 6 < Bq < 0. The last inequality follows by monotonicity of fi^ in 
and being a nonincreasing function of uj. For any w""" — 1 < z < y and 

0>0, 



' /(x) ^ fjx + l) ^ fiy + l) y-z+i 
-^i/fx + f) 11 e« - V e** / 

X—Z ^ ' X—Z 



Given < y < uj™^^ and 1 > e > 0, there is a large enough which makes the 
last fraction smaller than e. With such a choice we have 



l-e 



P'{^<y}= E /W</(y + i) E £'~'+'<^' 



Therefore, for the case of a finite oj"^^^^ choosing y — uj"^^^ — 1 and large 
makes (|B.3p arbitrarily close to oj"^^^. When oj™'*'' = oo, the argument shows 
that g{0) > y + I can be achieved for any y > 0. A similar computation 
demonstrates that any density towards w™'" can be reached when = —oo. □ 

Corollary B.3. The measures fi^ are stochastically nondecreasing in g. 

Proof. Since g and are strictly increasing functions of each other, it is equiv- 
alent to show monotonicity of fi^ . This follows if we can show < -^E^(^(w)) 
for an arbitrary bounded nondecreasing function ip. Lemma IB . II transforms this 
derivative into the covariance of (p{uj) and oj, which is non-negative due to ip 
being nondecreasing. □ 

Monotonicity of /i^ requires somewhat more of a convexity argument. 

Proposition B.4. The family of measures jl^, defined in (|2.17p . is stochasti- 
cally nondecreasing in g. 

Proof. Start by rewriting the definition: 

- g] ■ > y}) _ Cov^(w, 1{lu > y}) 



Vare(cj) Cov'^iu, uj) 



d 



P'{to > y} 



^-;Lpn.>y}. 
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Let us denote the /I^-expectation by E^. Fix a bounded nondecreasing function 
Lp. We need to show 

< A.Y,e^(^^Y 
ag 

We compute a different expression for this derivative. Passing the derivative 
through the sum in the third equahty below is justified because the series in- 
volved are dominated by certain geometric series, uniformly over 9 in small open 
neighborhoods. This follows from the definitions of 9_ and 6 and the assumption 

e_ < eig) < e. 

= E ^(2/)-^[P^{t^>y}-i{0>y}] 
= ^ E ^(2/) • [P^^ > 2/} - i{0 > y}] 

= -j-g^' E ^(y) • [li'^ >y}^ > y}] 

^Tg^' E f{y)-m^>y>o}~Ho>y>u;}] 



Above we introduced the function 

x-l 

<i>{x) = J2^iy) - E^(y)' 

y=l y=x 

with the convention that empty sums are zero. To conclude the proof, notice 
that $(a; + 1) — = f{x). Thus a nondecreasing function ip determines a 
(non-strictly) convex function <f> with $(1) ™ 0, and vice-versa. Hence Theorem 
lA. II establishes that 



_d 

dg ^ ' 



— E'?(^(w) = — E^$(w) > 0. □ 
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C Regularity properties of the hydrodynamic 
flux function 

For the zero range process defined among the examples in Section 12.21 the 
hydrodynamic (macroscopic) flux function H : K.+ — > M+ of (|2.14p is given by 



The results of Section [X] for / now read as follows: 

Proposition C.l. // the jump rate f of the zero range process is convex (or 
concave), then the flux Ti is also convex (or concave, respectively). Moreover, 
in this case Ti."{g) > (or Ti."{g) < 0, respectively) for all g > if and only if 
f is not a linear function. 

Parts of this proposition were proved with coupling methods in . 

Next we show in the general case that Ti.{g) is well defined, and is infinitely 
differentiable. (We use third derivatives in the proof of Theorem 12.31 ) The 
function Ti.{g) is, in general, the expected net growth rate w.r.t. fi^ as defined in 
(|2.14p . We show that the series making up this expectation is finite, even after 
differentiating its terms. This will then lead to smoothness oiTi{g). 

Lemma C.2. Letg{y, z) > be any function onZxZ, bounded by a polynomial 
in \y\ and \z\. Then for any 9 < 9 < 6, 

'E'^[{p{ujo, uji) + q{ujo, uji))g{ujo, LOi)] < oo. 

Proof. We deal with the first part that contains p, the one with q can be treated 
analogously. The sum we are looking at is 



These sums are certainly convergent if w""" and w™^'^ are both finite. When this 
is not the case we split both summations at zero, and convergence is established 
on the four quadrants of the plane. We use (|2.7p and the corollary 



n{g) = ESfiu;). 



UJ' 




f{y)\-f{z)\ • Z{9y 



fJ>(y+z) ]_ 



+ 1 z—uj- 




<y <uj 



,max 



and LJ' 



,min 



< Z < UJ 



,max 



p{y, z) < p{y, 0) = p{l, y-l)-^< p{l, 0) • 





and the corresponding part of the summation is bounded by 
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y < 0, ^; > 0: In this case 

p{y, z) < p{l, 0), 
and the corresponding part of the summation is bounded by 



• y < 0, z < 0- In this case 

p{y, z) < z) = p{z + 1,0)- -ffl- < p{l, 0) ^^^^ 



and the corresponding part of the summation is bounded by 



e{v+z) I 



p{l, 0)/(l) . J2 J2 g{y, z) . • 



2/ > 0, z < 0: In this case 



,(.,.)=P(. + 1,.-1)-;^<P(1,0).-^, 
and the corresponding part of the summation is bounded by 



P(1>0)- E E 9{y,z) 



Convergence of each of these bounds for ^ < ^ < ^ is estabhshed e.g. by the 
ratio test. □ 

Notice that a similar argument gives finite higher moments of the rates when 
log(/) is at most Unear in both directions on Z. 

Corollary C.3. H{q) is infinitely differentiable at all g e (u™", w™"^). 
Proof. By the previous lemma the series 

, max 

m ■■ = nm) = ^ • J^Jpiy^ ^) - i(y^ '^^jwiW- 

is convergent and infinitely differentiable. Since 7i(gi) = F(6{q)) and q i— > 6{q) 
is infinitely differentiable as well, the claim follows. □ 
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